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J>~>' Abstract 

I ' In this work, the hnear stabihty of the viscous incompressible fluid flow 

^ . between two parallel horizontal porous stationary plates with the assumption 

M-H ' that there is a small constant suction at upper plate and a small constant injec- 

C/3 . tion at the lower plate is studied. The Navier-Stokes and continuous equations 

are reduced to an equation modified by the suction Reynolds number, which 
^ . we call modified Orr-Sommerfeld equation. This equation is rewritten as an 

^H ' eigenvalue problem and is solved numerically using Matlab (Windows Ver- 

O.. sion). The effect of small suction Reynolds number on the linear stability 

fluid flow is discussed. 
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^ . modified Orr-Sommerfeld equation. 
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f^ ■ 1 Introduction 

^' 

O 
en 



The flow through porous boundaries is of great importance both in technologi- 
cal as well as biophysical flows. Examples of this are found in soil mechanics, 
the aviation industry, transpiration cooling, food preservation, cosmetic industry, 
blood flow and artificial dialysis. A large number of theoretical investigations deal- 
ing with steady incompressible laminar flow with either injection or suction at the 
C^ . boundaiies (||T1 -||2]) have appeared during the last few decades. Several authors 

||3l have studied the steady laminar flow of an incompressible viscous fluid in a 
two- dimensional channel with parallel porous walls. Rioual and al. (1996) stud- 
ied the power balance of a flat plate used as an airplane wing and found that there 
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is an optimum speed of suction that reduces drag and weaken energy consump- 
tion. Recent experimental studies show that suction damps growth of disturbances 
induced by free stream turbulence and transition is delayed/prevented (Fransson 
and Alfredsson 2003). Aspiration may also be used as a tool to induce transition, 
instead of the delay, if for example the wall material chosen is not able to provide 
continuous suction (MacManus and Eaton, 2000) or if it is applied non-uniformly 
(Roberts andFloryan, 2001) H. 

In this work, the flow of the incompressible viscous fluid between two paral- 
lel horizontal stationary porous plates with the assumption that there is constant 
small suction at upper plate and a constant small injection at the lower plate is 
considered. Speeds suction and injection are assumed uniform and same standard 
V* (see figure ([T])) below |i3|. The objective is to study the effect of number Reuj 
of extradvective term introduced in an On-Sommeifeld equation by small suction 
velocity on the stability of the fluid flow. Such attempt has been made earlier by 
E. Niklas Davidsson and L. Hakan Gustavsson the 2000s ( UJ and |4|) the case of 
the boundary layers (a flat plate-law) with unifomi wall suction by normalizing all 
the components of the velocity with the free stream velocity U*. Before E. Niklas 
Davidsson and al. ; V. M. Soundalgekar, V. G. Divekar and al. (1968 - 1972) stud- 
ied the flow quilt with suction at the stationary plate and laminar slip-flow through 
uniform circular pipe with small suction ( |I1]-||6]). All these studies confirm the 
stabilizing effect of small wall suction on a flow. 

The paper is organized as follows. In the second section the modified Orr-Sommerfeld 
equation governing the stability analysis in the Poiseuille two parallel horizontal 
porous stationary plates flow will be checked. In the third section an analysis of 
the small suction Reynolds number R^c^ stabilizing effect will be investigated. The 
conclusions will be presented in the final section. 

2 modified Orr-Sommerfeld equation 

The motion of a fluid is governed by the conservation laws, such as conservation 
of mass and momentum. The momentum of the fluid changes due to the forces 
acting upon it, and is governed by Newton's second law of motion. Fluids at low 
velocity compared to the speed of sound can be considered incompressible, i.e. the 
density of the fluid is constant. The stress acting on a fluid element may be linearly 
proportional to the local pressure and velocity gradient. For an incompressible, 
newtonian fluid the momentum and conservation equations yield 



dx* 



(2) 



when written by Einsteins' summation convention. These are the Navier-Stokes 
(N— S) equations. Here x* and u* are the i*^™^ component of the space and velocity 



vectors, as commonly in the literature denoted {x*, y* ,z*) and {u*,v*,w* 
streamwise, wall-normal and spanwise directions, respectively. Also, P* is the 
pressure, p is the density and u is the kinematical viscosity of the fluid. 

We Considered the viscous and incompressible fluid flow between two parallel 
horizontal stationary porous flat plates under a constant gradient of pressure with 
small wall suction and injection in the figure ^ below ||3l. 
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Figure 1: Poiseuille porous parallel flat plates flow. 

The coordinates (streamwise x, the wall-normal y and spanwise z) are scaled 
with the length scale /i* ( distance between the two walls). The streamwise and 
spanwise velocities u and w, respectively, are scaled with the streamwise free 
stream velocity U* while the wall-normal velocity v is scaled with the chai^acter- 
istic velocity of suction and injection V* . The pressure p is scaled with p[/*^ and 
the time t with jp. The Reynolds' numbers used here are defined as R^ = ^ ^ 



and suction Reynolds' number i?e^ = " where v is kinematic viscosity. 
We want to study the linear stability of a high Reynolds i?e and small suction 
Renynolds iJeoj (0 < Reu) < 1) numbers flow. The non-dimensional Navier-Stokes 
and continuous (Conservation of mass) linearized equations ([T]) and ^ where the 
variables are normalized as above take the following forms 



dt dx Re dy dz dx Re 

dv d^ Reuj_^dv_ dv _ __^^, J_ 2~ ... 

dt dx Re dy dz Reuj dy Re 

dw _ dw Reu) _ dw _ dui dp 1 2 ~ 

di dx Re dy dz dz Re 



du Rew_dv_ , ^ _ Q .gs 

dx Re dy dz 

For stability analysis, the flow is decomposed into the mean flow and the distur- 
bance according to 

Ui{r,t) = Ui{r)+Ui{r,t), (7) 

p{r,t) = P{r)+p{r,t). (8) 



We take the dimensional base flow for small suction and injection (see IS) 

,2 



base flow 



u{y) = u*{i-^) 


(9) 


V = v: 


(10) 


W = 0. 


(11) 


s above, we obtain with h* - 


= ±1 the no-dimensional 


U{y) = (l-y2) 


(12) 


V = 1 


(13) 


W = Q. 


(14) 



To obtain the stability equations for the spatial evolution of three-dimensional, we 
take the dependent on time disturbances 

(■u(x, y, z, t); v{x, y, z, t); w{x, y, z, t);p{x, y, z, t)); (15) 

which are scaled in the same way as above. 

By using continuity, the pressure terms can be eliminated from Navier-Stokes equa- 
tions. For such a mean profile (base flow), the divergence of Navier-Stokes equa- 
tions and continuity gives 

The equations © and (IT6l) after linearization give 



d 2 TT (^ 2 ^eui d 2 d?U dv _ 1 4 

dt dx Re dy dy"^ dx Re 

The disturbances are taken to be periodic in the streamwise, spanwise direc- 
tions and time, which allow us to assume solutions of the form 

f{x,y,Z,t) = /(y)e^("-+/3.-a.t). (ig) 

where / represents either one of the disturbances u, v, w or p and / the amplitude 
function, k, a = k^ = k cos 9 and /3 = ky = k sin 9 are the wave numbers, 
uo = ac the pulsation of the wave. With i"^ = —1, 9 = {kx,k), c = Cr + ici 
wave velocity which is taken to be complex, a and /3 are real because of temporel 
stability analysis considered. 
Then with the equation (fTSl ). the equation ([17] ) becomes 

ia[{U-c){D^-k^)-U"]v = -^D{D^-k^)v 



+^ {D' - k')' v; (19) 



where D = -^ ; with boundary conditions for all (x, ±1, z, i) 

v{±l) = 1 
v'{±l) = 0. (20) 

Taking 

vp{x,y,z,t) = v{y)e'(--+P^~-^)-l, (21) 

the equation (fT9l ) and the boundary conditions take the forms 

(u - ia-i^D^ (D^ - k") - U" - ia-'^{D' - k^f] Vp = 

\ tie J Re J 

c (L»2 - k^) Vp, (22) 



^p(±l) = 

v'p{±l) = 0. (23) 

The equation (l22l) is a flow equation modified by suction Reynolds number Re^j 
(or the speed of suction and injection), which we call modified Orr-Sommerfeld 
equation, rewritten as an eigenvalue problem. 



3 Stability analysis 

We consider our three-dimensional disturbances. We use a temporal stability anal- 
ysis as mentioned above. With c complex as we have defined above, when Cj < 
we have stability, q = we have neutral stability and else we have instability. We 
employ Matlab (Windows Version) in all our numerical computations to find the 
eigenvalues. The Poiseuille parallel horizontal stationary porous plates flow with 
the basic profile 

U=(l-y2,l,0) (24) 

for Reuj small ( i.e. small suction) is considered (see lO). The eigenvalue problem 
(l22l) is solved numerically with a suitable boundary conditions. The solutions are 
found in a layer bounded at y = ±1 with U(ibl) = (0,1,0). The results of 
calculations are presented in the figures below. We present the figures related to 
the eigenvalue problem (l22l ). 

For all these figures the black, red, green and blue colors are respectively, for 
Reuj = 0, Reu} = 0.5, Reuj = 0.75, i?ew = 1 and the yellow color is for q = 0. 

In each group of four figures, the first one is Figure a), the second is Figure b), 
the third is Figure c) and the fourth is Figure d). 

For a fixed /c = 1, we get figure ([3]l of Cj vs. R^ for sequential values of 6. 
figure a) for 6 = 0, figure b) for 6 = O.lvr, figure c) for 9 = 0.2tt and figure d) for 
e = O.Svr. 



For a fixed k = 1.02, we get figure dUl of ci vs. R^ for sequential values of 9. 
figure a) for ^ = 0, figure 6) for 6 = O.lvr, figure c) for 9 = 0.2tt and figure d) for 
9 = O.Svr. 

For a fixed A: = 2, we get figure © of ci vs. i?e for sequential values of 9. 
figure a) for ^ = 0, figure 6) for 9 = O.lvr, figure c) for 9 = 0.2it and figure d) for 
= 0.37r. 

For a fixed A; = 3, we get figure ^ of ci vs. Re for sequential values of 9. 
figure a) for 9 = 0, figure 6) for = O.Itt, figure c) for 9 = 0.2tt and figure d) for 
61 = O.Svr. 

Through the Figures (l3]and|4]i, it is easy to see that the stability increases when 
Reui increases for all 9 because, for any familly of curves of these figures the slope 
fall when the suction Reynolds' number increases. The Reynolds' critical number 
Rec for which we have the transition becomes important when i?e^ increases (see 
tabular Q) which confirms that the wall small suction or injection have a stabi- 
lizing effect on the viscous incompressible flow. Because of R^ui influences, we 
said that to normalize with the characteristic velocity of suction is necessary for 
a perfect command of the field of stability. In particular for (Reu = 0,k = 1.02, 
without suction, figure dU black curve (a) ) we find i?ec = 5772 which corresponds 
exactly to the critical value given by classical linear theory for a plane-Poiseuille 
flow without suction or injection. 

For i?e^ = 1, the equation (l22l) is identicaly for the equation (2.28) found by 
E. Niklas Davidsson and L. Hakan Gustavsson in HI. For this value of R^c^ in the 
cases of figures (l3]and|4ll the field of stability is more greater than the other small 
suction Reynolds number for all 9. 

Note that for great k the case of the Figures ^ and Q there is stability without 
transition for all value of 9 and Reco but when the small suction R^c^ increases the 
field of stability decreases. We can also say that k influences the stability of the 
flow i.e. the growth of the wave number induces also the stability of the flow. We 
therefore conclude that in two parallel horizontal stationary porous plates viscous 
and incompressible fluid flow, with small suction and small injection, the small 
suction Reynolds number stabilize the flow for small wave number but the high 
wave number effect is important than small suction Reynolds' number on the fluid 
flow stability. 
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Figure 2: Tabular for critical Reynolds numbers 
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Figure 3: Growth rate q vs. Reynolds' number for A: = 1. 
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Figure 4: Growth rate q vs. Reynolds' number for k = 1.02. 
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Figure 5: Growth rate q vs. Reynolds' number for k = 2. 
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Figure 6: Growth rate q vs. Reynolds' number for k = 3. 

4 Conclusion 

In this paper, we have investigated the effect of small suction Reynolds number 
on the stability of the fluid flow between two pai^allel horizontal stationary porous 
plates . We have shown that the instability of the perturbed flow is governed by a 
remarkably equation named modified Orr-Sommerfeld equation. We find also that 
the normalization of the wall-normal velocity with characteristic small suction (or 
small injection) velocity is important for a perfect command of fluid flow stability 
analysis. We noticed that for k = 1.02 and Reuj = we find Rec = 5772 which 
corresponds exactly to the critical value given by classical linear theory for a plane- 
Poiseuille flow without suction or injection. We noticed also that the high wave 
number stabilize more than the small suction Reynolds number. 
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